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Abstract. We have surveyed a sample of 165 solar-type spectroscopic binaries (SB) with periods from 1 to 30 days 
for higher-order multiplicity. A subsample of 62 targets were observed with the NACO adaptive optics system 
and 13 new physical tertiary companions were detected. An additional 12 new wide companions (5 still tentative) 
were found using the 2MASS all-sky survey. The binaries belong to 161 stellar systems; of these 64 are triple, 11 
quadruple and 7 quintuple. After correction for incompleteness, the fraction of SBs with additional companions 
is found to be 63% ± 5%. We find that this fraction is a strong function of the SB period P, reaching 96% for 
P < 3 d and dropping to 34% for P > 12 d . Period distributions of SBs with and without tertiaries are significantly 
different, but their mass ratio distributions are identical. The statistical data on the multiplicity of close SBs 
presented in this paper indicates that the periods and mass ratios of SBs were established very early, but the 
periods of SB systems with triples were further shortened by angular momentum exchange with companions. 
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1. Introduction 

The formation of binary stars remains a subject of active 
research and debate (Zinnecker & Mathicu 2001). Close 
binaries are particularly difficult to explain because at or- 
bital periods of a few days, the separations are less than 
the size of the individual components during the proto- 
stellar phase. Fission of a rapidly rotating pre-main se- 
quence stellar configuration appears unlikely according to 
the most recent theoretical results (Tohline 2002). Thus, 
some mechanism(s) of orbit shrinkage must be present in 
order to extract the angular momentum from a proto- 
binary. This momentum could be deposited in wider stel- 
lar companions - a hypothesis we test in this paper. 

Stars form in groups and interact dynamically. 
Detailed studies of dynamical decay of small unstable 
clusters, e.g. by Sterzik & Durisen l|1998f) . were success- 
ful in explaining such properties as multiplicity rate and 
mass ratio distribution. But dynamical interactions within 
star clusters alone cannot sufficiently broaden an initially 
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narrow separation distribution (Kroupa & Burkert l2001|) . 
Additional dissipative processes and even higher stellar 
densities are required (Sterzik, Durisen & Zinnecker l2004fl . 

Larson 1|20021) argues that tidal torques excited by stel- 
lar companions in an accretion disk are the most likely 
and efficient mechanism to extract angular momentum 
from accreting stars. In the absence of angular momen- 
tum transport, accreting matter accumulates in a disk 
or torus that fragments rapidly, creating the companion. 
Indeed, hydrodynamical simulations of Bate et al. ( 2002 ) 
show how such companions are formed and how accreting 
binaries become tighter by interacting with their distant 
companions or with fly-by members of a nascent cluster. 
Reipurth ( 2000 ) gives convincing evidence that strong ac- 
cretion and jet activity is actually observed in embed- 
ded young multiple systems. He relates periodic knots 
in Herbig-Haro jets to periastron passages in inner (un- 
resolved) binaries and interprets the decreasing separa- 
tion between the knots as a fossil record of orbit shrink- 
age. However, the triple systems examined by Reipurth 
and Bate et al. are relatively wide systems (~ 10 A.U.), 
and cannot account for the characteristics found in short 
period pre-main-sequence (PMS) spectroscopic binaries 
(SB). Several PMS binaries with periods shorter than 2 
days are known (Melo et al. [2001), an d the SB fraction 
among PMS stars may be as high as among field stars. 
Many PMS SBs are actually members of higher-order mul- 
tiple systems (Sterzik et al. I2004fl . A comprehensive star 
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formation theory must be able to explain these observa- 
tional facts. 

The discovery of massive extrasolar planets with pre- 
dominantly short orbital periods revived the interest in 
orbital decay via interactions within disks. Planet migra- 
tion in a disk is now a generally accepted theory. A distant 
(stellar) companion enhances planet growth and migration 
(Zucker & Mazeh EES Udry et al. 12171131 Eggenberger et 
al. |2)04 ) , possibly in an analogous way to multiple-star 
formation as envisioned by Larson H2002[) and others. 

Present-day parameters of close binaries are not iden- 
tical to their parameters at birth (i.e. when their masses 
were assembled) because of subsequent evolution. One 
such evolutionary mechanism is called KCTF - Kozai 
cycles with tidal friction (Eggleton & Kisseleva-Eggleton 
2001 ). A distant companion in a stable (hierarchical) triple 
system causes periodic modulation of the inner-binary ec- 
centricity by Kozai cycles. These cycles are periodic, but 
only as long as the inner system does not interact tidally 
at periastron. In this case the Kozai cycles are gradu- 
ally modified, the inner system gets "locked" in a high- 
eccentricity state and its orbit then slowly decays to a 
circular SB with a period of few days (Kisseleva et al. 
1998 ). The period distribution of solar-type close binaries 
within higher-order multiple systems seems to match this 
scenario, showing a sharp drop in the number of systems 
with P >7 d (Tokovinin & Smekhov l2lJU2) . 

A competing evolutionary effect is the disruption of 
multiple systems by dynamical instabilities or by pertur- 
bations from other stars passing nearby. Thus, SBs formed 
within higher order multiples can also lose tertiary com- 
panions (TCs) at later times. Interactions with field stars 
can be neglected in the present study because they disrupt 
only very wide tertiaries with separations above 0.1 pc 
(see Close et al. 1199 for a review of observational data 
and theory). Much closer tertiaries with periods as short 
as ~ 300 yr could be "ionizied" during early evolutionary 
stages in very young clusters, explaining the difference in 
wide-binary frequency between young T Tau associations 
and the field (KroupaEUJUl)- 

Observers have noticed the high multiplicity of close 
binaries for some time. Mayor & Mazeh l|1987fl looked for 
orbital precession (presumably caused by TCs) in a sam- 
ple of 25 solar-type SBs and found that some 25% of those 
binaries are triple. Isobe et al. (1992) searched for visual 
companions to SBs by means of speckle interferometry. 
Tokovinin l|1997l MSC) showed that 43% of nearby solar- 
type stars with periods under 10 days have known ter- 
tiaries. All five such systems in the Duquennoy & Mayor 
(|1991l DM91) G-dwarf sample are triple. 

Are all close binaries triple? Is the angular momentum 
of the close binary system always extracted by a tertiary 
component (via the KCTF mechanism or during accre- 
tion)? In the present paper, we attempt an observational 
investigation of the question and prove that not all close 
SBs are triple. 

The paper is organized as follows. We describe our 
sample of close solar- type binaries in Sect. [21 New compo- 



nents were discovered with adaptive optics (Sect. EJ) and 
the 2MASS sky survey (Sect.^J. We estimate the massses 
and periods of all tertiaries in Sect. and describe the 
detection limits of new and existing techniques in Sect. [SJ 
We study the statistics of tertiary components in Sect. [7J 
Section |H1 summarizes the results, and the implications for 
close- binary formation theories are discussed in Sect. El 

2. The sample 








B-V 

Fig. 1. Absolute V- magnitudes of the SB targets in the 
main sample versus their B — V colors. Asterisks mark 
the SBs with evolved components, the solid line shows 
the Main Sequence from Lang (1992). 

We define a sample of spectroscopic binaries (SBs) 
where the chance of discovering a TC is maximized. 
Ideally, the systems should be composed of low- mass dwarf 
stars because such stars are not too bright, increasing the 
possibility of detection of close visual companions, and 
have sharp lines making it easier to detect spectroscopic 
tertiaries. Furthermore, the stars should be nearby, lead- 
ing to larger separation of tertiaries. We also want to probe 
the frequency of TCs to SBs with periods below and above 
the KCTF limit of P - 7 d , and hence consider the SBs 
with periods up to 30 d . 

The main sample was constructed by selecting solar- 
type SBs from the catalog of Batten et al. {1989). No 
condition on the luminosity class was imposed: short or- 
bital periods guarantee small stellar radii, hence evolved 
components (giants) will be excluded. It is known that 
luminosity classes assigned by observers are often contra- 
dictory, hence we tried to ignore them. It turned out later 
that some objects are in fact moderately evolved (Fig. [TJ. 

The sample was complemented by new SBs from recent 
surveys of nearby stars (Latham et al. 120021 Goldberg et 
al. 2002, Halbwachs et al. 2003|and other). New unpub- 
lished SBs from the CORALIE ongoing survey were also 
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Table 1. Preliminary parameters of CORALIE binaries. 
The columns contain: (1) HIP number, (2) HD number, 
(3) orbital period, (4) semi-amplitude of the primary and 
(5) Notes. 



HIP 


HD 


Pi, 


Ki, 


Note 






d 


km/s 




2790 


3277 


15 


5? 


Uncertain SB 


2848 


3359 


22 


22.4 


B-comp. to HIP 2888 


19248 


26354 


2.5 


56 


P 3 > 3000 d 


38179 


64184 


17.8 


12.7 




45957 


81044 


15 


8 




49161 


87007 


30 


13.4 


P3 suspected 


53217 


94340 


6.8 


20? 


Susp. P 3 = 1200 d 


56960 


101472 


4.4 


15.7 




73269 


132173 


14 


31.5 




85675 


158577 


9.69 


34.5 




94863 


180445 


2.5 


47 


SB2 


97030 


186160 


10.7 


18 




107779 


207450 


4.9 


15.8 




114703 


219175 


7.1 


22.9 


CPM HIP 114702 


116429 


221818 


11.6 


15 





added (Table [TJ. A total of 200 SBs were thus selected. 
However, the main sample was then restricted to stars 
with Hipparcos parallaxes larger than 10 mas (i.e. within 
100 pc). This sample contains 165 SBs with periods less 
than 30 days belonging to 161 stellar systems (including 
four multiple systems with two short-period SBs each). 
Each system is identified by its Hipparcos number (when 
several components have distinct HIP numbers, only one 
is retained) . The main data on the sample - HIP and HD 
numbers, coordinates, proper motions, parallaxes, magni- 
tudes and spectral types - are given in Table 2, available 
electronically. Several targets belong to the Hyades clus- 
ter. 

The masses of the primary components of the SBs are 
cither known from orbital solutions (e.g. eclipsing pairs) 
or estimated from the B — V colors and spectral types us- 
ing standard relations for the Main Sequence (MS) from 
Lang (1992). Despite the brightness of our targets, the 
task of assembling the basic data turned out to be non- 
trivial because the light of several components is mixed 
and there is a confusion between components in the cat- 
alogs. For 32 stars identified as evolved, we roughly es- 
timated their masses Mi from the absolute magnitudes 
V abs as Mi = 0.2(9 - V abs ). Those stars in Fig. [2 that 
are above the MS but not marked with asterisks belong 
to multiple systems that contain more massive evolved 
visual components, while the spectroscopic primaries are 
unevolved. The median mass of a primary in our sample 
is Mi = I.IMq, the full range is 0.42 to I.IMq. 

The masses of spectroscopic secondaries M2 are calcu- 
lated from the known mass ratios in case of double-lined 
systems (SB2s) or from the mass functions in case of SBls 
(minimum masses). It is known that minimum masses are 
statistically biased (Goldberg et al. 2003). However, we 



are not concerned with the mass ratio of SBs, so this bias 
does not affect the present investigation.. 

3. Adaptive-optics observations 

Adaptive optics (AO) imaging in the near infrared (IR) is 
a powerful tool for discovering low-mass companions, as 
demonstrated by a number of recent studies (e.g. Shatsky 
& Tokovinin 2002 ) . AO allows for an exploration of a large 
part of the parameter space that is not accessible to tra- 
ditional observational techniques. This increased capac- 
ity results from the high angular resolution, high dynamic 
range (the ability to detect faint stars because AO concen- 
trates light into diffraction-limited cores) and observations 
in the infrared which improves the contrast of low-mass 
companions with respect to their primary stars. 



6669 


11843 


19248 




• • 




21395 


25760 


41834 



Fig. 2. Representative narrow-band NACO images of 
some resolved targets in 2" x 2" fields (inverted intensity 
scale with square- root stretch). The Hiparcos numbers are 
marked. 

High spatial resolution images of the target stars have 
been obtained at the VLT telescope using the NAOS- 
Conica adaptive optics system, NACO 1 , on November 8-9, 
2004 and on July 6-12, 2005. We observed 72 targets and 
two astrometric calibrators (see below) in the narrow-band 
(NB) filter centered at the wavelength 2.12 /im, to avoid 
detector saturation. Five stars were also observed in the 
J and/or H photometric bands to estimate the colors of 
their companions. Only targets without known TCs (or 
with distant tertiaries) were observed. Of the 72 stars ob- 
served, only 62 belong to our final sample. The remaining 
ten SBs (HIP 5436, 25760, 34003, 35600, 41834, 50966, 
87330, 87370, 110514, 112009) do not satisfy the criteria 
for inclusion in the study sample, but we nevertheless pro- 
vide data on their companions in the electronic table. 

The images were processed in a standard way. The sky 
background and detector bias were estimated and sub- 
tracted by median filtering a series of five frames dithered 
on the sky. The de-biased images were then flat-fielded and 

1 http:/ /www.eso.org/instruments/naco 
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combined. The relative positions and magnitudes of wide 
components were determined by DAOPHOT point spread 
function (PSF) fitting of the main star to all sources. 
Measurements on individual images indicate that the ran- 
dom error in position is about 0.5 mas in each coordinates 
and the rms error in magnitude difference is around 0.02 m 
for components with Am < 3 m , increasing to 5 mas and 
0.05 m for Am = 5 m . 

Data on close (< 1") companions were processed by 
an alternative technique of simultaneous fitting of binary 
parameters and a non-negative PSF. First, we determine 
the initial binary parameters in Fourier space by fitting 
(at high spatial frequencies above the seeing cutoff) the 
image Fourier transform (FT) to the FT of a PSF star con- 
volved with the binary. This initial estimate of the binary 
parameters then serves to de-convolve the object FT from 
the binary, extracting an estimate of the true (synthetic) 
PSF. The synthetic PSF is radially-averaged outside the 
ten pixel (0'.'13) radius and is used in the second, final 
model-fitting. In this way, it is possible to measure binary 
parameters independently of the PSF star. The quality of 
this procedure is controlled by the abscence of companion 
traces in the synthetic PSF before its radial averaging. 

Most of the new companions are well above the de- 
tection threshold (Fig. |2J). In trying to recover faint and 
close tertiaries that escape immediate detection, we exam- 
ined all images after subtracting radially-averaged profiles. 
This procedure reveals asymmetric features of the PSF, of- 
ten persistent in the individual (non-averaged) images. We 
measured the parameters of eight such tentative compan- 
ions and found that their separations are all close to O'.'l, 
with a magnitude difference of about 2 m — 3 m . The posi- 
tion angles of these companions differ by 170° to 180° from 
the VLT parallactic angle, strongly suggesting that they 
are telescope-related artefacts. In the case of HIP 107095, 
a very good radial- velocity coverage (Feke HTM7|) virtually 
excludes the existence of our tentative O'.'l TC, unless its 
orbital inclination is very low. In contrast, the new close 
companion to HIP 19248 is real: it does not follow the 
parallactic-angle dependence and is confirmed by other 
techniques. 

Relative coordinates of the components in detector 
pixels were transformed to on-sky positions using ob- 
servations of two known wide binaries, HIP 108797 and 
HIP 116737. The accuracy of this calibration is entirely 
determined by the known angular separation and posi- 
tion angle of these "calibrators" . We adopted the result- 
ing pixel scale of 13.30 ± 0.01 mas/pixel and offset of 
—0.4° ± 0.1° to the observed angles. The calibration re- 
mained stable between 2004 and 2005 runs. The data on 
newly discovered tertiaries are presented in Table 01 We 
also include the data on our calibrator binaries in the last 
rows. We thus detected 13 new TCs (as well as the known 
component HIP 98578B) among 62 targets from the main 
sample (detection rate 21%). For some new TCs, we were 
able to check the correspondence of their colors with those 
expected for the MS dwarfs (photometric confirmation), 
while two TCs are confirmed by proper motions. We es- 



Table 3. Data on companions discovered with NACO. 
The columns give (1) Hipparcos number, (2) companion 
identification, (3) spectral band, (4) epoch of observation, 
(5,6,7) measured separation, position angle and magni- 
tude difference. 



HIP 


Cmp 


Band 


Epoch 


P, 


0, 


Am 








2000+ 


II 


O 




6669 


B 


NB 


4.857 


1.469 


69.37 


5.50 






J 


4.857 


1.467 


70.31 


6.29 






H 


4.857 


1.466 


70.04 


5.84 


11843 


B 


NB 


4.857 


1.422 


339.2 


3.55 


19248 


B 


NB 


4.857 


0.104 


46.11 


2.61 






H 


4.857 


0.094 


47.18 


2.39 


21395 


B 


NB 


4.857 


0.991 


5.15 


3.61 


25760 


B 


NB 


4.857 


0.751 


184.20 


4.21 


41834 


B 


NB 


4.857 


1.047 


337.93 


3.85 






J 


4.857 


1.056 


336.47 


4.51 




C 


NB 


4.857 


3.573 


202.83 


5.30 


43557 


B 


NB 


4.857 


3.446 


333.73 


4.28 


44164 


B 


NB 


4.857 


1.851 


270.83 


6.93 






J 


4.857 


1.855 


270.80 


7.29 


48215 


B 


NB 


4.857 


0.759 


124.09 


2.49 


64219 


B 


NB 


5.533 


0.309 


291.41 


4.05 


91360 


Opt 


NB 


5.526 


3.760 


31.76 


7.59 


94863 


B 


NB 


4.857 


9.378 


52.20 


4.68 


98578 


C 


NB 


5.517 


0.391 


340.23 


2.25 


98578 


B 


NB 


5.517 


3.696 


353.64 


0.62 


107779 


B 


NB 


4.857 


2.219 


144.60 


3.56 






J 


4.857 


2.226 


144.59 


4.16 






NB 


5.522 


2.233 


144.30 


3.45 






Calibrators 








108797 


B 


NB 


4.857 


3.848 


246.37 


2.52 


116737 


B 


NB 


4.857 


3.890 


276.96 


0.98 



timate that all new companions except two are physical, 
and support this with statistical agruments presented in 
Sect. 4. Comments on selected systems are provided be- 
low. 

HIP 19248 is an astrometric binary in the Hipparcos 
catalog. We estimate the TC period as ^6 yr, its 
mass as 0.3 Mq- The TC was previously inferred from 
CORALIE's precision RVs. 

HIP 35487 has a known astrometric TC with com- 
puted orbit (P3 = 92 yr, estimated separation 0'.'4, mass 
0.34 A4q) which was not detected with NACO; it must 
be a white dwarf. 

HIP 41834 has two companions in the NACO images. 
However, the fainter companion C was not detected in the 
J-band, thus it must be an optical source. 

HIP 91360 has such a faint companion that, if physi- 
cal, it would have sub-stellar mass. We consider this com- 
paion an optical source because the density of background 
objects near this target is very high, TV* = 335 (Sect.^J). 

HIP 94863: a wide companion was seen by luck in the 
corner of a dithered field. It has been independently iden- 
tified as a 2MASS photometric candidate (Sect.^J, with a 
consistent position. The proper motion would have moved 
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it by Of! 7 in 5 yrs if it were only an optical source, hence 
we conclude that it is physical. 

HIP 98578 is a visual triple, with the new AC pair dis- 
covered here. We estimate the period of AC to be ~100 yr; 
it could be detectable as astrometric perturbation in the 
motion of AB. The 2356-yr visual orbit computed for AB 
('Hopmann ll973[l is premature. Without AO observations, 
we would wrongly consider the B component as tertiary. 

HIP 107779. The new companion is confirmed as phys- 
ical by our second-epoch observation and its J — K color. 
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4. Search for companions in 2MASS 
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Fig. 3. CMDs for HIP 93966 in a crowded (N* = 187) 
field (top, two optical candidates) and HIP 66514 with 
A* = 8 and two physical companions (bottom). The pri- 
mary targets are marked by large squares. The selected 
photometric candidates are plotted as large asterisks, re- 
maining field stars as crosses. Stars within 30" from the 
target are marked by small squares. The solid line depicts 
the standard Main Sequence from Lang (1992). 



The SBs were examined for the presence of wide vi- 
sual companions using the Two Micron All-Sky Survey 
(2MASS) and the Digital Sky Survey (DSS). The primary 
targets are bright nearby dwarfs. Many faint "compan- 



100 



Separation, arcsec 



Fig. 4. Distribution of photometric candidates from 
2MASS in separation p and magnitude J. Known physi- 
cal companions are marked by squares, the new physical 
companions by diamonds, while crosses are likely optical 
(i.e. not physical) companions. 



ions" are found near any bright star, but the vast major- 
ity of such companions are simple line-of-sight projections 
or "optical" sources. However, faint physical companions 
may still hide among the multitude of optical compan- 
ions. Previous generations of double-star observers were 
guided by statistics to separate likely physical companions 
from chance projections (Poveda et al. ll982[) . Thus, a "fil- 
ter" was adopted to reject optical companions, and such 
a filter introduces a strong selection effect in the exist- 
ing double-star catalog, WDS (Mason et al. 2001). In this 
study, we use the JHK S photometry provided by 2MASS 
to place the companion stars on the color-magnitude dia- 
gram (CMD) and to select potential physical companions. 
The extinction in the J and K bands can be safely ne- 
glected for our nearby sample. 

The photometry and coordinates of all stars within 
a 2' radius from each of the 165 targets were retrieved 
from the 2MASS Point Source Catalog using the Vizier 
service at CDS. A total of 6079 companions were thus 
found. Only 202 (3.3%) were selected automatically as 
photometric candidates by their distance from the MS 
y/A J 2 + A( J - K) 2 < 0.2 m . We checked that this crite- 
rion selects most of the known physical companions, but it 
cannot be made "sharper" without the risk of losing com- 
panions. Two examples of CMDs are shown in Fig. [3J It 
is clear that in the fields with a large number of stars A* 
the confusion impedes the detection of physical compan- 
ions and that the photometric candidates in these fields 
are likely optical. 

The statistics of photometric candidates (Fig. con- 
firms that they are mostly optical - faint and with sep- 
aration p of the order of the field radius, 120". In con- 
trast, known physical companions are brighter and are 
distributed in logp almost uniformly. 

In an effort to reduce the number of candidates, we 
discarded 36 fields with A* > 40 (the median is N* = 13). 



6 



Tokovinin et al.: Tertiary companions to SB 



Faint companions with p > 60" and J > 13 m were also 
discarded, as well as multiple candidates and candidates in 
the Hyades (HIP 20019, 20284, 20440). Excluding known 
physical companions, we retained 35 candidates for astro- 
metric checking. 

The targets are nearby stars with typical proper mo- 
tions (PMs) of O'.'l per year. Physical companions must 
share these motions and would be displaced by ~ 5" on 
the sky in 50 yr. Old photographic plates are over-exposed 
in the vicinity of the targets, hence we can check only can- 
didates with p > 20". The images of fields around each 
target were retrieved from the DSS (red Palomar plates 
for northern declinations and red UK Schmidt plates for 
southern declinations). The time base between 2MASS 
and DSS is around 50 yrs in the North but only around 
10 yrs for southern stars. 

We found photometric candidates in the DSS images 
and roughly measured their coordinates using the astro- 
metric reference of the DSS, with an accuracy of about 
0"5. When the displacement of the main target is signif- 
icantly larger than 0"5, we can discriminate the optical 
companions. The majority of wide and faint companions 
turned out to be optical. In four cases the status of the 
companions remained uncertain (small PM and/or small 
time base), while three new companions with large sepa- 
rations (HIP 16042, 36328, 66514) among our 35 selected 
candidates were confirmed. 

The distribution of new companions in the (log p, J) 
plane (Fig. HJ shows a "cluster" of 5 points near p 10". 
These TCs have not been checked by astrometry. However, 
statistical arguments show that they are all likely physi- 
cal companions. We leave out HIP 94863 which lies in a 
crowded field (iV* = 52), but note that its 9" 4 compan- 
ion is physical (Sect. The remaining 4 stars with such 
companions are HIP 3362, 17076, 20712, 60331. These tar- 
gets have small stellar surface densities (iV* of 17, 13, 15, 
6, average 12.75) and the companion separations are less 
than 11". Among 86 targets with N* < 13, the expected 
number of point sources at p < 11" is 8.7. Given that the 
photometric criterion selects only w3% of all sources, the 
expected number of randomly selected photometric candi- 
dates within 11" is 0.26; this number should be compared 
to the four companions actually found. In addition, there 
are 6 previously known physical TCs with similar separa- 
tions of p ~ 10". 

An additional screening for new close companions has 
been performed by examining the 2MASS images of ap- 
parently single targets. We found such companions to 
HIP 43557 (confirmed with NACO) and HIP 60956. The 
latter is still considered as tentative. 

Summarizing the results for 2MASS candidates, 
Table 01 lists the data on 23 known, 7 new certain, and 
5 new tentative (marked with ?) companions - a total of 
35 TCs. The known companions are marked "MSC" . We 
see that existing catalogs contain about 2/3 of the wide 
TCs; another 1/3 are added in this work. 



Table 4. Tertiary companions found in 2MASS. The 
columns give (1) HIP number, (2,3) separation and po- 
sition angle, (4) J magnitude of the companion, (5) J — K 
color, (6) status, with MSC for previously known compan- 
ions. 



HIP 


P," 


0,° 


J 


J — K 


Status 


3362 


11.18 


278.13 


11.366 


0.801 


new 


7874 


10.45 


31.34 


8.705 


0.749 


MSC 


12189 


38.11 


274.44 


6.080 


0.258 


MSC 


16042 


95.82 


129.08 


8.182 


0.556 


new,CPM 


17076 


9.89 


245.34 


11.857 


0.839 


new 


20712 


11.03 


119.67 


11.795 


0.861 


new 


24663 


9.95 


148.60 


8.481 


0.460 


MSC 


31850 


30.17 


248.74 


11.420 


0.840 


MSC 


36238 


43.00 


323.59 


8.249 


0.780 


new,CPM 


41211 


35.49 


328.69 


11.441 


0.841 


new? 


42172 


10.20 


25.11 


5.882 


0.385 


MSC 


45957 


8.25 


118.02 


10.179 


0.787 


MSC 


47053 


24.95 


149.20 


7.142 


0.271 


MSC 


52064 


3.67 


342.14 


6.136 


0.034 


MSC 


56809 


9.06 


248.40 


6.573 


0.615 


MSC 


56960 


79.81 


334.28 


10.403 


0.705 


MSC 


60331 


10.19 


320.43 


12.921 


0.896 


new 


60956 


5: 


355: 


- 


- 


new? 


61910 


58.52 


228.80 


8.822 


0.444 


MSC 


66514 


56.34 


258.63 


7.662 


0.532 


MSC 


66514 


65.86 


254.25 


11.610 


0.882 


new,CPM 


73269 


81.52 


39.64 


9.800 


0.716 


new? 


74037 


17.80 


19.75 


9.239 


0.494 


MSC 


84586 


33.97 


92.78 


8.542 


0.913 


MSC 


91009 


16.91 


45.50 


10.424 


0.902 


MSC 


94863 


9.48 


53.03 


10.457 


0.769 


NACO 


98578 


3.00 


353.70 


6.251 


0.103 


MSC 


103569 


10.55 


67.95 


6.419 


0.283 


MSC 


104026 


39.99 


201.07 


12.735 


0.865 


new? 


107354 


14.56 


289.05 


8.224 


0.572 


MSC 


108461 


13.69 


246.71 


5.984 


0.379 


MSC 


111802 


24.52 


350.36 


7.344 


0.853 


MSC 


114379 


15.65 


75.94 


6.002 


0.332 


MSC 


114639 


87.20 


205.54 


11.033 


0.805 


new? 


114703 


25.12 


355.67 


6.484 


0.338 


MSC 



5. Parameters of tertiary companions 

Table \5\ contains the data on previously known compan- 
ions to the SBs which have been extracted from the MSC 
(Tokovinin 1997) and combined with the new detections 
in a single database. Additional bibliographic searches for 
stars without reported tertiaries have been made through 
SIMBAD. We also scanned the fourth catalog of interfer- 
ometric measurements (Hartkopf et al. 2004 - INT4) and 
marked the stars that were observed by speckle interfer- 
omctry. Also, the ninth catalog of spectroscopic binary 
orbits (Pourbaix et al. 120041 - SB9) was searched for new 
orbits. 

A part of Table [S] with the information on all known 
components is printed here; the full Table is available elec- 
tronically. Below we describe the information in this Table. 
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Table 5. Data on companions (fragment). The columns are explained in the text. 



HIP 


ID 


Type 


Separat 


ion 


Period 




Mi 


M 2 


Remark 


47053 


Aa,Ab,A 


S2 


0.471 


mas 


1.681 


d 


1.54 a 


1.46 a 


SB 575 


47053 


A,C,AC 


V 


0.300 




64.257 


y 


3.00 s 


0.98 a 


COU 2084 


47053 


AC,B,* 


Chrp 


25.000 




41.099 


ky 


3.98 s 


1.65 s 


STF 1369 


47053 


Ba,Bb,B 


SI 


2.530 


mas 


28.231 


d 


1.29 a 


0.36 m 


1998AstL...24..288T 


48215 


Aa,Ab,A 


SI 


0.924 


mas 


3.390 


d 


1.11 a 


0.29 m 


SB 584 DI Leo 


48215 


A,B,* 


AO 


0.759 




130.282 


y 


1.40 s 


0.61 v 


NACO-2004 


48273 


Aa,Ab,* 


S2 


1.196 


mas 


3.055 


d 


1.21 a 


1.15 q 


SB 585 


48833 


Aa,Ab,* 


S2 


1.057 


mas 


3.055 


d 


1.26 a 


1.20 q 


SB 584 


49018 


Aa,Ab,A 


E,S2 


0.690 


mas 


1.070 


d 


0.78 a 


0.53 q 


DH Leo, SB 591 


49018 


A,B,* 


V 


0.220 




14.190 


y 


1.31 s 


0.50 v 


CHARA 145 


49161 


Aa,Ab,* 


SI 


3.965 


mas 


30.000 


d 


0.74 a 


0.50 q 


CORALIE 


49809 


Aa,Ab,* 


SI 


5.557 


mas 


28.098 


d 


1.49 a 


0.20 m 


SB 601 



HIP 47053 



AC 



Aa 



Aa.Ab.A P=1.68d 



Ab 



Ba 



Bb 



A,C,AC P~64yr 



AC,B,* P~41000yr 



Ba.Bta.B P=28.2d 



Fig. 5. Example of the hierarchical quintuple system 
HIP 47053 and the identifications of its sub-systems (see 
text). 



Additional comments on specific systems are provided in 
the electronic notes to Table [5] 

Identification of each system is given in the column 2 as 
a sequence of 3 designations (primary, secondary, parent) 
separated by commas. These codes describe the hierarchy 
of each system by referencing to the parent, as illustrated 
in Fig. [SJ The system at the highest hierarchy level (root) 
is coded with an asterisk, it is the parent of the wide sub- 
system AC,B in Fig. |SJ 

The type code of each system shows the method(s) of 
its discovery, in the same manner as in the MSC. For ex- 
ample, SI and S2 refer to single- and double-lined SBs, 
E to eclipsing systems, C stands for wide (CPM) sys- 
tems (the following small letters describe which criteria 
of physical relation between components are satisfied, see 
MSC and notes to the electronic Table 5), v means a re- 
solved system closer than 3", V stands for a visual binary 
with computed orbit, etc. Two new special types are AO 
(companions discovered with NACO) and 2M (2MASS 
companions). 

The separation between components is given together 
with units (arcseconds or milliarcseconds). For systems 
with known orbits the separation refers to semi-major 
axis, otherwise it is the observed separation p. The separa- 
tion of SBs is estimated from known periods and Kepler's 
Third Law as 



P 



(1) 



where p is the parallax (arcseconds), P is the orbital pe- 
riod in years and the system mass A4 is in solar units. 
Here we assume implicitly that separation is statistically 
equivalent to semi-major axis. 

Periods are given in the next column, in units of days, 
years or kiloyears. For wide systems without computed 
orbits the periods are estimated from p with Eq. 

Masses of primary and secondary components Ai\ and 
M.i are given in solar mass units with codes indicating 
the method of mass estimation. The preferred methods 
are a - from spectral type (cf. Sect. |2J or * - estimated 
by other authors from orbit solutions or detailed mod- 
els. Other codes are 1 - evolved components, s - sum of 
sub-components' masses, q - SB2 secondary, m - mini- 
mum mass of SB1 secondary, v - mass from magnitude 
difference, : - unknown magnitude difference, ? - uncer- 
tain companions. Masses of all new components (AO and 
2M types) were estimated from the magnitude differences 
in the K band using the relations of Henry & McCarthy 
(1993) and taking into account the light of the spectro- 
scopic secondary. 

Remarks provide common identifications for the sub- 
systems, e.g. the SB numbers in Batten et al. (0-9893 cat- 
alog, visual double-star designations, bibcodes, etc. 

The sample contains 165 spectroscopic binaries belong- 
ing to 161 stellar systems; 79 have no known tertiaries, 
and the remaining 86 SBs have one or more companions 
(64 triples, 11 quadruples, seven quintuples - a total of 
82 systems). Of these, six TCs are uncertain (our 2MASS 
detections, doubtful speckle companions, etc.), but they 
are still considered in the statistics. The fraction of SBs 
with TCs should be no less than 86/165=0.52. Not all 
TCs have been discovered yet. Known distant compan- 
ions are considered as tertiaries, but in some systems the 
true tertiaries at closer separations are not known and the 
present "tertiaries" are actually higher-level companions, 
as turned out to be the case for HIP 98578, where a new, 
closer companion to SB was discovered with NACO in 
addition to the known wide companion. 
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Mass ratios of TCs are defined here as q% — 
i.e. relative to the spectroscopic primary, not to the to- 
tal mass of the SB. There is no ambiguity for triples, 
but for systems with more components, as in Fig. |3J the 
choice of tertiary is not always evident. We select M3 
to be the largest mass of all components at the next 
higher level of hierarchy. Thus, the binary Ba,Bb in Fig. [3] 
has Aa as its most massive tertiary, (73 = 1.10/1.29 = 
0.85, P3 « 41000 yr. However, another close binary 
Aa,Ab in the same system has a tertiary component C, 
hence q 3 = 0.98/1.10 = 0.89, P 3 » 64 y. 

6. Detection limits 



1 .0 H 



0.8 - 



0.6 



0.4 



0.2 - 



0.0 




Fig. 6. Estimated limits of companion detection by differ- 
ent techniques for HIP 64219. Full line - radial velocities, 
dotted line - NACO, dashed line - visual and CPM, dash- 
dotted line - 2MASS. The NACO companion is marked 
by asterisk. The overall probability of TC detection for 
this object is zero in the hatched zone below all curves 
and one above. 



For any given star and any observing technique, the 
probability of detecting a TC depends on the compan- 
ion's parameters - its period, P3, and its mass ratio, q 3 . 
This probability should be a smooth function of the pa- 
rameters, but we simplify it here as a sharp limit, assum- 
ing that all TCs with q 3 > qu m (P3) are detected and less 
massive companions are not. An example of such limits 
for one target is shown in Fig. |S| By averaging the detec- 
tion probabilities for the whole sample, we again obtain 
a smooth detection probability (Fig. 0). The method of 
computing detection limits is detailed in Appendix A for 
each observing technique. We specify the detection lim- 
its for each system, hence the four systems with two SBs 
each have only one entry per system. However, these sys- 
tems are at least quadruple and pose no problem in the 
statistics. 

Detection limits were converted into estimates of the 
minimum observable mass ratio (73 (Sect.|3J and compared 
to the actual mass ratios of known and new TCs in Fig. EI 




Fig. 7. Probability of tertiary companion detection for the 
whole sample as a function of log P3 and (73 . 



100 
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10 
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"0 b 



Fig. 8. Comparison of the actual mass ratios of tertiary 
components 93 with their estimated detection limit qu m . 
Known TCs are plotted as pluses, the new (NACO and 
2MASS) TCs as asterisks. Uncertain TCs are marked with 
squares. 

We see that only several tertiaries are slightly below the 
estimated detection limit and conclude that the limits can 
be used confidently in the statistical analysis. 

7. Statistics of tertiary companions 

7.1. Mass ratios and periods of tertiaries 

Figure shows the distribution of TCs in ^3,(73) space. 
Known companions at higher hierarchical levels (like the 
component Ba relative to the Aa, Ab system in Fig. |5j) 
are also plotted as squares for comparison. The curves 
indicate detection limits. As expected, new TCs detected 
in this work mostly have low mass ratios. 

We note that 16 tertiaries (17 ± 4%) have (73 > 1, i.e. 
the spectroscopic primary is not the most massive compo- 
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Fig. 9. Mass ratio of tertiary to spectroscopic primary as 
a function of tertiary orbital period. The known compo- 
nents are plotted as pluses, the new components as aster- 
isks. The squares denote components at the next hierar- 
chical level, not considered in the main statistical analysis. 
The lines trace detection probabilities of 10% (lower) and 
50%(upper), i.e. the contours of the surface in Fig.0 



nent in the system. If components were selected randomly, 
one would expect this fraction to be 1 /3 or higher (in sys- 
tems with more than 3 components we select the most 
massive tertiary). Hence, we confirm the known tendency 
that the most massive components in multiple systems are 
preferentially found in close sub-systems. 

Interestingly, we find only three TCs with J3 > 
10 5 yrs, but also six higher-level companions wich such 
long periods. Thus, a decrease of the number of TCs at 
long periods must be real, not a selection effect. 

7.2. SBs with and without tertiary companions 

The cumulative distributions of the orbital periods of SBs 
with and without TCs are compared in Fig. QJJI It is ev- 
ident that SBs within multiple systems have, generally, 
shorter orbital periods. The maximum difference between 
the cumulative distributions normalized to 1 reaches 0.30. 
The numbers of objects in the histograms are 86 and 
79, so the Kolmogorov-Smirnov test rejects the hypoth- 
esis that both distributions are equal with a significance 
level of 0.999. Some of the SBs considered today as "sin- 
gle" have yet undiscovered tertiaries (e.g. suspected as- 
trometric companions), so the actual difference between 
the distributions is likely to be larger. The presence of a 
TC shortens the period of a spectroscopic binary. On the 
other hand, the existence of such difference shows that 
SBs without TCs (pure binaries) do exist. 

In contrast, there is no difference in the mass ratio dis- 
tributions of the SBs with and without TCs: both show 
a slight excess of systems with nearly equal-mass compo- 
nents (twins), but otherwise the distributions are almost 
uniform (Fig. 1110. As noted in Sect. [21 the mass ratios of 




Period, d 

Fig. 10. Histograms and cumulative distributions of or- 
bital periods of spectroscopic binaries with tertiaries 
(Triple, full line) and without tertiaries (Binary, dashed 
line). The larger fraction of short periods among those 
SBs which belong to triples is evident. 




Fig. 11. Cumulative distributions of the mass ratios of 
spectroscopic binaries with tertiaries (Triple, full line) and 
without tertiaries (Binary, dashed line). 

SB1 are biased, but this bias is independent of the pres- 
ence of the TC. 

7.3. Period-period diagram 

There seems to be no correlation between the SB or- 
bital periods Pi and the TC periods P3 (Fig. aster- 
isks). However, the shortest TC periods (lower bound- 
ary of points) seem to increase with P\. For comparison, 
we overplot nearby solar-type multiples from the MSC 
(crosses, cf. Tokovinin 2004). Multiple systems with wide 
inner binaries are almost exclusively found with period 
ratios P3/P1 of between 5 and 10 000. In contrast, our 
multiples are all quite far from the dynamical stability 
limit and can have period ratios > 10 4 . This difference 
may be understood if the SBs were formed with longer pe- 
riods and later migrated to shorter periods. The relative 
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P,, d 

Fig. 12. Comparison of the SB periods (horizontal axis) 
with respective tertiary periods (vertical axis). The full 
line denotes the dynamical stability limit P3 = 5-Pi, the 
dotted line marks the period ratio of 10 000. Our sample 
is plotted as asterisks, other triple systems from MSC are 
plotted as crosses. 

paucity of inner sub-systems with periods of a few months 
matches this scenario, although it could probably also be 
explained by observational selection. The lowest observed 
P3 ~ 2 yr implies that the progenitors of these SBs had 
Pi < -P3/5 ~ 150 d , otherwise those systems would have 
disintegrated through dynamical instability (e.g. Mardling 
& AarsethUlTJ). 

7.4. Companion distribution by Maximum Likelihood 



two roughly equal parts: close (Pi < 7 d, 90 systems) and 
wide (Pi > 7 d, 75 systems), and repeat the ML analysis 
for these sub-samples. The results are shown in Table HO 
and in Fig. The errors of the TC frequency are deter- 
mined by the shape of the S-function near its minimum, 
as illustrated in Fig. IB.l| We also provide in Table El the 
raw (uncorrected) TC frequencies / raw to show that our 
correction for incompleteness is not dramatic. The calcu- 
lation was repeated ignoring uncertain TCs, and only a 
slight reduction of the estimated TC frequency results. 

0.25 [ 1 1 1 1 1 



P,<7d 
P,>7d 




Pj. Y r 

Fig. 13. Period distributions of TCs around close SBs 
(squares) and wide SBs (triangles). The error bars corre- 
spond to ±1(7. The distribution of solar-type binaries from 
DM91 is plotted for comparison (dotted line and crosses). 



An estimate of the TC distribution over period P3 is 
derived here by the Maximum Likelihood (ML) method. 
The parameter space is divided into bins, the fraction of 
TCs in each bin is f = {/&}. The probability of obtaining 
the actual data C is called the likelihood. By minimizing 
its natural logarithm S 

S(f) = -21n£(f), (2) 

we find the estimate of f. Moreover, the ML method 
permits us to determine confidence intervals of the esti- 
mated parameters or their functions (Avni et al. I198()|l . 
Mathematical details and ML equations are given in 
Appendix B. 

In the analysis, we select five bins for P3 uniformly 
covering the logP3 space from one to 10 6 yrs. The de- 
tection probabilities are computed as a surface above the 
<7iim(P3) curve in each period bin. This is strictly true only 
for a uniform distribution of (73. The actual 53 distribution 
does in fact appear rather uniform (Fig.[5J). A more correct 
analysis of the joint TC distribution in P3, 53 (with 3 bins 
in 53) has been performed and gives practically identical 
results. 

The total TC fraction estimated by ML, i.e. corrected 
for incompleteness, is 63 ± 5%. Considering the difference 
between short- and long-period SBs, we split our sample in 



Table 6. Frequency of tertiary companions in two sub- 
samples. N is the number of systems in each sample, / raw 
is the TC frequency uncorrected for incomplete detections, 
/ is the TC frequency estimated by the ML method. The 
numbers in italics are obtained by discarding uncertain 
companions. 



Sample 


N 


/raw 


/ 


Pi < 7d 


90 


0.66 


0.80 ± 0.06 




88 


0.65 


0.79 ± 0.06 


Pi > 7d 


75 


0.33 


0.40 ± 0.06 




71 


0.30 


0.36 ± 0.06 


All 


165 


0.52 


0.63 ± 0.05 



The sample was then further sub-divided into four 
groups ranked over SB periods, and the calculation of 
the TC frequency was done for each group separately 
(Fig. [Til). For the shortest periods, Pi < 2.9 d , the esti- 
mated TC frequency is as high as 96%. 

It is clear that many SBs have tertiary companions. 
But companions are typical also for normal solar-type field 
stars (Duquennoy & Mavor lT991[l . Is there any difference 
between SBs and field stars with respect to wide compan- 
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Fig. 14. Frequency of tertiary companions as a function 
of SB period. 

ions? We compare in Fig.EDthe companion frequencies in 
these samples, both corrected for selection effects. Overall, 
the companion frequency in the DM91 sample for periods 
above one yr is 50% - higher than in the wide-SB sub- 
sample. On the other hand, the fraction of companions to 
close SBs is clearly increased. Their periods range from a 
few years to several thousand years. 

7.5. High multiplicities 

Our sample contains 64 triples and many systems of multi- 
plicity higher than 3: 11 quadruples and seven quintuples 
(Sect. 5). Batten l|1973fl introduced /„ - the fraction of 
systems of multiplicity n and higher to systems of mul- 
tiplicity n-l. For our sample, fa = (H + 7)/64 = 28% 
and fa = 7/11 = 64%. These numbers are comparable 
to the high- multiplicity frequencies in the MSC: fa = 
179/626 = 29% and fa = 38/141 = 27% (Table 1 of 
Tokovinin 12001(1 . Spectroscopic binaries from the MSC, 
such as HIP 76563/76566 (two SBs of 3.27 d and 14.28 d 
in a system with six components), are not included in the 
sample. This is done to avoid intentional bias towards high 
multiplicity. 

8. Summary 

1. The period distributions of SBs with and without 
tertiary companions (TCs) are significantly different. 
This proves the existence of pure SBs without TCs. 

2. The mass ratio distributions of SBs with and without 
TCs are identical. 

3. The frequency of TCs is 63% for the whole SB sam- 
ple. However, it is a strong function of the SB period, 
reaching 96% for the close (Pi < 3 d ) SBs and decreas- 
ing to 34% for SBs with P x > 12 d (Fig. Hi). The last 
number is less than the companion frequency of ^50% 
for G-dwarfs in the field. These results are robust be- 
cause the TC detection does not depend on the SB 
period. 



4. The periods of most TCs in our sample range from 2 yr 
to 10 5 yr. There is no correlation between P3 and Pi. 
The triple systems with Pi < 30 d have large period 
ratios and are very stable dynamically fFig. 1120. 

5. The TCs are more massive than spectroscopic pri- 
maries in a small fraction (17% ± 4% ) of multiple 
systems, suggesting that there is a tendency of the 
most massive component to be found preferentially in 
shortest-period sub-systems. 



9. Discussion 

It seems that the properties of close SBs are established 
early in their evolution. In a compilation of pre-Main 
Sequence (PMS) SBs by Melo et al. ijMITj) . the period 
distribution is not different in any significant way from 
the distribution in our sample. For example, SB periods 
are divided almost equally between l d — 7 d and 7 d — 30 d 
intervals. The observed frequency of TCs in the Melo et 
al. PMS sample, 42% ± 19% (Sterzik el at. EUMjl . is not 
different from the observed TC frequency of 52% in our 
sample of SBs in the field. 

Our main result (Fig. I14fl shows that essentially all 
very close SBs are members of higher-order systems. The 
KCTF mechanism is most likely responsible for shortening 
SB periods during their life on the Main Sequence. It acts 
independently of the TC period and mass: even a very 
wide and low-mass companion can influence the SB, given 
sufficient time. The period of Kozai cycles is of the order 
of P^Pz/P-i) ~ 10 6 yr for a typical tertiary period of 
P3 = 1000 yr and the period of an SB progenitor of Pi = 
1 yr. It takes many (> 10 2 ) cycles to complete the KCTF 
evolution, hence it could not occur at the PMS stage. The 
formation of very close PMS SBs is not possible because 
the PMS stellar radii are too large and a contact PMS 
binary would rapidly merge (Whelan 0.970). 

For SB periods longer than 3 days, we find an increas- 
ing proportion of pure SBs without tertiaries. The exis- 
tence of non-triple SBs is clearly demonstrated by this 
study. We see that the frequency of wide TCs to bina- 
ries with Pi > \2 d is lower than the companion frequency 
to single solar-type stars. These pure SBs have not passed 
through KCTF evolution. What, then, was the mechanism 
which produced their orbital decay? 

Sterzik, Durisen & Zinnecker 123224} discuss various 
processes of binary formation in a systematic way and 
provide suitable scaling laws. Molecular cores fragment 
into mini-clusters of the size of order 100 AU. The system 
crossing time turns out to be shorter than the accretion 
time scale, hence stars build their masses and interact dy- 
namically simultaneously. An accreting multiple system 
becomes tighter (Umbreit et al. I2005fl . Close approaches 
of the components of an unstable multiple system lead 
to ejections and formation of a close binary with a semi- 
major axis ~ 10 times shorter than the system size. Joint 
action of accretion and disruption can thus produce close 
SBs without additional companions. To form a 10 d SB 
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(semi-major axis 0.1 AU), we need a typical precursor 
multiple system of ~1 AU size. 

Not all primordial multiples are unstable. If orbits of 
distant companions contain the bulk of a system's angular 
momentum (as should often be the case), there are no close 
approaches causing disintegration. The hardening of inner 
orbits in such systems can be primarily driven by accretion 
and interactions with a disk. An interesting alternative has 
been suggested by Namouni (2005): an orbit can become 
eccentric if the system moves with acceleration, e.g. caused 
by an asymmetric jet or a TC, and then tidal forces will 
produce a close SB in analogy to the KCTF process. 

Some stable TCs will be lost at later stages by dy- 
namical interaction with other members or systems of a 
primordial star cluster (Kroupa 2001), again leaving pure 
SBs. 

Accretion on a close binary biases its mass ratio to- 
wards 1. Indeed, the mass-ratio distribution of SBs is 
nearly uniform (Goldberg et al. [2003 ) . SBs with equal- 
mass components (twins) can be explained this way 
(Tokovinin EUnOI Bate & Bonnell UMTjl . It is remarkable 
that the mass ratio distributions of SBs with and without 
TCs are the same (Fig. lll|) . pointing to a common process 
of their mass buildup. A general tendency to find the most 
massive component of a multiple system as a primary in 
the SB is confirmed here. It can be explained both by ac- 
cretion physics and by pure iV-body dynamics. However, 
we find that distant TCs have often comparable masses or 
are even more massive than the SB (Fig. |5J), contradict- 
ing the simulations of Delgado-Donate et al. (J2003J) who 
predict only low-mass TCs. 

It becomes thus increasingly clear that close binaries 
were formed by a combination of different processes (some 
probably not yet identified) and that their orbits and 
multiplicity were further modified by subsequent evolu- 
tionary effects such as KCTF. How efficient is the KCTF 
mechanism in producing close SBs? Case-by-case model- 
ing (Kisseleva et al. I1998fl should now be complemented 
by a statistical analysis. Additional information on the 
KCTF will be obtained from the relative orientation of 
orbits in close triple systems which are now observable 
with long-baseline interferometry (e.g. Muterspraugh et 
aL EDDSjl . 

This work can be considered as a small step towards 
complete multiplicity statistics in the solar neighborhood, 
still largely unknown because of strong observational se- 
lection effects. For example, is the paucity of ~ 100 d inner 
periods in multiple stars (Fig. I12|) real? Is it a fossil record 
of inner-orbit decay? A systematic census of all compan- 
ions in nearby solar-type stars using a combination of ob- 
serving techniques is required to answer these questions. 
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Appendix A: Discussion of detection limits 

The knowledge of the detection limits of various observing 
techniques is essential to characterize the completeness of 
the TC statistics (Sect. 6) and to correct the final results 
(Sect. 7.4). The limits of each technique are discussed and 
modeled below. 

A.l. Adaptive optics 

We studied the detection limits for triple companions in 
our NACO observations of spectroscopic binaries. At each 
radial distance from the spatially unresolved SB the rms 
fluctuation a in the flux over the corresponding circle was 
computed, and the detection limit was assumed to be 5a. 
This procedure has been checked by simulating artificial 
companions and then detecting them visually. The result- 
ing detection limit is described by a model 

AK < 7r, p< O'/l 

AK < 5r + 0.9, O'/l < p < l'/5 

AK < 9.05, p>l"5, (A.l) 

where r — log(p/0'/035). The diffraction limit at A = 
2.12 fim is X/D = 0'.'054. A similar model was used by 
Shatsky & Tokovinin l|2002|l . Our NACO detection limits 
are corroborated by Brandeker l|20041 his Fig. 4). 

A. 2. Visual and CPM companions 

The detectable magnitude differences of known visual 
companions increases with their separation p. The upper 
envelope corresponds to 

AU < 4.51og(p/0'/05). (A.2) 

The components lying at this limit were all discovered 
with speckle interferometry (cf. INT4). Thus, (|A. 2|) may 
be somewhat optimistic for 1/2 of our systems that have 
never been observed with speckle techniques according to 
INT4. No new companions are given in the CHARM cat- 
alog of lunar occultaions and long-baseline interferometry 
by Richichi et al. ( 2004). An additional constraint that the 
apparent magnitude of a companion should be brighter 
than U = 15™ is introduced because fainter stars are 
generally not detected as common-proper-motion (CPM) 
companions. 

A. 3. Radial velocities 

A TC modulates the center-of-mass velocity of the SB 
with a longer period. The probability of detecting such 
modulation depends on the precision of radial-velocity 
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observations a, their number N and the time span AT. 
Studies of detection limits with simulations (Halbwachs et 
al. 2003) established that the relevant parameter is the ra- 
dial velocity amplitude Aq of a circular-orbit binary with 
the same period and 90° inclination, 



213M 3 Pn 



-1/3 



(A.3) 



where Aq is in km/s, tertiary mass M3 and the SB system 
mass A4 are in solar masses, and the tertiary's period 
P3 is in days. Following Halbwachs et al., we model the 
spectroscopic detection limit conservatively as A /a > 3 
for all periods shorter than the data span AT, and as 
A /a > 3 + 201og(T 3 /AT) for periods shorter than 2 AT. 
The actual detection limit may extend to periods P3 much 
longer than AT if the data are analyzed carefully, but we 
prefer to be conservative here by assuming the sharp loss 
of detection capacity at P3 > 2AT. 

For each SB, we determined the relevant parameters 
a and AT either from the data given in SB9 (Pourbaix 
et al. I2UU4|) or from the original publications. The quality 
of spectroscopic orbits varies greatly, from crude orbital 
solutions made some 80 years ago (e.g. HIP 5689, 72524, 
86263) to high-precision velocities (HIP 80686) or long- 
term coverage (HIP 107095). Moreover, orbits are often 
computed by combining radial velocities of different qual- 
ity or with different systematic offsets. Thus, any model 
of spectroscopic detection limit is necessarily very crude 
and Aq/<j is only an indicative parameter, at best. 

Despite uncertainties in the spectroscopic detection 
limits, we note that half of our sample has AT > 5.5 yr, for 
78% of SBs AT > 2 yr. The median error is a = 1 km/s. 
The TCs with periods of a few years would have a good 
chance of being discovered spectroscopically, hence the 
paucity of such companions in our sample must be real. 

A. 4. Astrometry 

Astrometry is a poweful technique to detect low-mass 
companions with orbital periods over several years. 
Makarov & Kaplan (2005.) recently published a catalog 
of astrometric binaries revealed either by a significant 
difference between "instantaneous" proper motion (PM) 
measured by Hipparcos and long-term PM in the Tycho-2 
catalog (A/i binaries) or by the acceleration measured di- 
rectly by Hipparcos (/i binaries, also known as G-solutions 
in Hipparcos). There are 11 astrometric binaries in com- 
mon with our sample, of which six are already listed here 
with TCs discovered by other techniques (including HIP 
19248 and 64219 resolved for the first time here). We did 
not resolve two ji binaries, HIP 35487 and 114639, with 
NACO, possibly because their astrometric companions are 
white dwarfs (cf. discussion in Notes to Table 5). 

Gontcharov et al. I|2001[> published a similar study: 
they brought historical astrometric catalogs into the 
Hipparcos system and identified several stars as astro- 
metric binaries. However, only a few astrometric orbits 



of these stars are published to date (none for our sample) , 
and Makarov & Kaplan ( 2005 ) did not independently con- 
firm the PM variations of many stars. We observed with 
NACO but did not resolve HIP 67153, 80686, 85365, 86263 
listed as binaries by Gontcharov et al. Here we consider 
astrometric detections without published orbits only as 
hints of the existence of TCs. Hence we do not take into 
account astrometric companions in our statistics. 

Appendix B: Maximum Likelihood method 

Problem layout. Suppose that a total of N targets have 
been surveyed. We determine the distribution of compan- 
ions over parameters by dividing the parameter space into 
K bins, the fractions of companions in each bin are de- 
noted as f = {fk}- Let Ni t k be the number of detected 
companions for i-th star and fc-th bin (1 when a compan- 
ion is detected, otherwise). The total number of com- 
panions per bin is Nk = J2i Nik- Thus, a raw estimate of 
companion frequency will be 



/fc,raw = N k /N. 



(B.l) 



In fact the probabilities of companion detection for 
each star and each bin d^k are generally less than 1, so 
we have to account for incomplete detections. This is best 
done by the Maximum Likelihood (ML) technique. First, 
we evaluate the probabilities of observations for each star 
Pi and compute the likelihood function 



(B.2) 



By minimizing the natural logarithm S = — 21n£ over 
parameters f we estimate the parameters. The contours of 
S in the f-space define confidence limits of the parameter 
estimates: AS — 1 corresponds to the 68% interval ("ler"), 
AS" = 2.71 to 90% and AS = 4 to 95%, in direct analogy 
with the Gaussian probability distribution (Avni et al. 

mm- 

Independent detections. In the simplest case we assume 
companion detections in each bin to be independent of 
other bins. The average expected number of detected com- 
panions is fkdi^k, and the probabilities are computed from 
the Poisson distribution, 



Pi 



(fkdiJ 



\N i>k 



■ exp {-fkdi,k}- 



(B.3) 



Considering that N^k only takes values of or 1, the like- 
lihood equation can be written as 

N K 

-21n£ = -2^^[^ fc ln(/ fc d a )-/ fc ^ fc ]+const.(B.4) 

i=l fe=l 

By differentiating (|B.4|I over fk we obtain the estimates 



fk 



B.5) 
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Fig. B.l. The shape of the S — S m i n function near its min- 
imum and the confidence intervals on the total companion 
frequency for a sub-sample of close SBs, Pi < 7d. 

This result is intuitive: we simply correct the observed 
companion frequencies / raw by the average probabilities 
of detection J2idi,k/N. 

We compute the probabilities di^ as fractional sur- 
faces of each bin above the detection limit curve qum(P3), 
individually for each star. The underlying assumption is 
that the (unknown) companion distribution within each 
bin is uniform. Looking at Fig. |3 we note that there are 
no strong gradients of companion frequency over or P3 , 
hence the assumption is good. 

Mutually exclusive detections. A given spectroscopic 
binary can have no more than one tertiary companion, 
by definition. Thus, a detection of a companion in sonic 
bin automatically excludes companions in other bins. 
Denoting by /o the probability of no companions, we have 
an additional condition fk + /o = 1. 

According to the Bayes theorem, the probability of ob- 
taining an observational result for the i-th star is related 
to the distribution fk and the conditional probabilities 
(detecting or not detecting companion) p(i\k), 

Pi = 5>(i|*)A. (B-6) 

fc=0 

In the case of companion non-detection, the probabil- 
ity p~ is 

Pi =fo+ E^ 1 - d ^f k = 1 - E ( BJ ) 

k=l k=l 

while for stars with companions the probabilities are sim- 
ply p\ — fkdi,k, according to (|B.3(I . This leads to the 
likelihood function 

-21n£ = -2J2H^--J2 d ^fk) 

i- k 

-2 E E N *,h H^.kfk) + const, (B.8) 

i k 



where i~ means that the summation extends over stars 
without detected companions. The second term reduces 
to a single sum over detected companions because other- 
wise Ni^ = 0. The parameters fk are estimated by direct 
minimization of l|B.8|) . These ML estimates are very close 
to those given by eq. IB. 51 

Confidence intervals. In order to find confidence inter- 
vals on fk or on the total companion frequency / = ^2 k fk, 
we have to perform constrained minimization of (|B.8|) with 
a fixed parameter /|. The curve S(f£) — S m i n then defines 
the confidence intervals, as illustrated in Fig. IB.l| 

References 

Avni, Y., Soltan, A., Tananbaum, L., & Zamoran, C. 1980, 
ApJ, 238, 800 

Bate, M.R, Bonnell, LA, & Brown, V. 2002, MNRAS 336, 
705 

Batten A.H. 1973, Binary and multiple systems of stars 

(Oxford: Pergamon). 
Batten, A.H, Fletcher, J.M, & McCarthy, D.G. 1989, Publ. 

Dom. Ap. Obs, 17, 1 
Bate, M.R. & Bonnell, LA. 1997, MNRAS, 282, 33 
Brandeker, A. 2004, in: Star Formation at High Angular 

Resolution, Proc. IAU Symp. 221, eds. M. Burton, R. 

Jayawardhana, T. Bourke, IAU, p. 229 
Close, L.M, Richer, H.B, & Crabtree, D.R. 1990, AJ, 100, 

1968 

Delgado-Donate, E.J, Clarke, C.J, & Bate, M.R. 2003, 

MNRAS, 342, 926 
Duquennoy, A. & Mayor, M. 1991, A&A, 248, 485 (DM91) 
Eggenberger, A, Udry S, & Mayor M. 2004, A&A, 417, 353 
Eggleton, P. & Kisseleva-Egggleton, L. 2001, ApJ, 562, 101 
Fekel, F.C. 1997, AJ, 114, 2747 

Goldberg, D, Mazeh, T, Latham, D.W. et al. 2002, AJ, 124, 
1132 

Goldberg, D, Mazeh, T, & Latham, D.W. 2003 ApJ, 591, 397 
Gontcharov, G. A, Andronova, A. A, Titov, O. A, & 

Kornilov, E. V. 2001, A&A, 365, 222 
Halbwachs, J.L, Mayor, M, Udry, S, & Arenou, F. 2003, 

A&A, 397, 159 

Hartkopf, W.I, Mason, B.D, Wycoff, G.L, & McAlister, H.A. 
2004, Fourth Catalog of Interferometric Measurements 
of Binary Stars, http://ad.usno.navy.mil/wds/int4.html 
(INT4) 

Henry, TD. & McCarthy, D.W. 1993, AJ, 106, 773 
Hopmann, J. 1973, Astron. Mitt. Wien, No. 13 
Isobe, S, Noguchi, M, Ohtsubo, J. et al. 1992, PNAOJ, 2, 459 
Kisseleva, L, Eggleton, P, & Mikkola, S. 1998, MNRAS 200, 
292 

Kroupa, P. 2001, in: Zinnecker, H, & Mathieu, R.D., eds. The 
formation of binary stars. Proc. IAU Symp. 200, IAU. 

Kroupa, P. & Burkert, A. 2001, ApJ, 555, 945 

Lang, K.R. 1992, Astrophysical data. Planets and Stars. Berlin: 
Springer- Verlag 

Larson, R.B. 2002, MNRAS, 332, 155 

Latham, D.W, Stefanik, R.P, Torres, G. et al. 2002, AJ, 124, 
1144 

Makarov, V.V. & Kaplan, G.H. 2005, AJ, 129, 2420 
Mardling, R.A. & Aarseth, S.J. 2001, MNRAS, 321, 398 
Mason, B.D, Wyckoff, G.L, Hartkopf, W.I. et al. 2001, AJ, 
122, 3466 



Tokovinin et al.: Tertiary companions to SB 



15 



Mayor, M. & Mazeh, T. 1987, A&A, 171, 157 
Melo, C.H.F., Covino, E., Alcala, J.M., & Torres, G. 2001, 
A&A, 378, 898. 

Muterspraugh, M.W., Lane, B.F., Konacki, M. et al. 2005, 

A&A, in print ( astro-ph/0509855 ) 
Namouni, F. 2005, AJ, 130, 280 

Nordstrom, B., Mayor, M., Andersen, J. et al. 2004, A&A, 419, 
989 80 

Pourbaix, D., Tokovinin, A. A., Batten, A.H. et al. 2004, A&A, 
424, 727 

Poveda, A., Allen, C, & Parrao, L. 1982, ApJ, 258, 589 
Reipurth, B. 2000, AJ, 120, 3177 

Richichi A., Percheron, I., & Khristoforova, M. 2004, A&A, 
431, 773 

Shatsky, N. & Tokovinin, A. 2002, A&A, 383, 92 
Sterzik, M.F. & Durisen, R.H. 1998, A&A, 339, 95 
Sterzik, M.F., Durisen, R.H., & Zinnecker, H. 2003, A&A, 411, 
91 

Sterzik, M.F., Melo, C.H.F., Tokovinin, A. A., & van der Bliek, 

N. 2005, A&A, 434, 671 
Tohline, J.E. 2002, ARA&A, 40, 349 
Tokovinin, A. 1997, A&AS, 124, 75 (MSC) 
Tokovinin, A. A. 2000, A&A, 360, 997 

Tokovinin A. 2001, in: The formation of binary stars. IAU 
Symp. 200, ed. H. Zinnecker & R.D. Mathieu, ASP Conf. 
Ser., 84 

Tokovinin, A. 2004, Rev. Mex. Astron. Astrof. Conf. Ser., 21, 
7 

Tokovinin, A. A. & Smekhov, M.G. 2002, A&A, 382, 118 
Umbreit, S., Burkert, A., Henning T. et al. 2005, ApJ, 623, 940 

2005ApJ...623..940U 
Udry, S., Mayor, M., & Santos, N. C. 2003, A&A 407, 369 
Whelan, J.A.J. 1970, MNRAS, 149, 167 

Zinnecker, H., & Mathieu, R.D., eds. The formation of binary 

stars. Proc. IAU Symp. 200, IAU. 
Zucker, S. & Mazeh, T. 2002, ApJ, 568, L113 



